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Non-Markovian master equation for interacting qubits coupled to a bosonic bath: 

analytic form and asymptotic approximation 
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Non-Markovian dynamics of two interacting two-level qubits coupled to a bosonic bath was pre¬ 
viously studied using the quantum-state-diffusion (QSD) equation, where a stochastic state is used 
to describe the system. In this study, we provide another perspective on this system by deriving the 
analytic form of the master equation, which describes the system with a reduced density matrix. 

Then, we validate the master equation by examining entanglement generation and state purity. In 
addition, with the master equation, we observe the effects from first-order noise and notice that a 
good asymptotic approximation to the master equation can be made by neglecting the first-order 
noise. 
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I. INTRODUCTION 

In the real world, a quantum system will always inter¬ 
act with the external environment and cannot be treated 
as an isolated system, leading to the introduction of 
open quantum systems [1]. Recently, open quantum sys¬ 
tems have attracted significant attention, and phenom¬ 
ena such as quantum decoherence [2-4] and quantum dis¬ 
sipation [5-7] have always been a key focus of several re¬ 
searchers in this field. Moreover, open quantum system 
concepts have been extensively applied in quantum op¬ 
tics [8], quantum measurement and metrology [3, 9-11], 
quantum control[12, 13], and quantum information sci¬ 
ence including quantum computation [14] and quantum 
simulation[15, 16]. 

Typically, an open quantum system is described by 
a reduced density matrix, and its dynamics is mathe¬ 
matically governed by a master equation. In a Markov 
approximation, the time evolution of a system is given 
by the well-known Lindblad master equation [17]. How¬ 
ever, the Markov limit is not always a good approxima¬ 
tion since open systems in Markov limit usually lose in¬ 
formation to the surroundings and exhibit irreversible 
dissipation and decoherence, which are bottlenecks to 
address in quantum information science. On the con¬ 
trary, non-Markovian process-based systems [1, 18, 19], 
where back reaction of environments is considered, re¬ 
gain the information previously lost to the surroundings. 
This makes the non-Markovian processes more physical 
and appropriate than the Markovian process under cer¬ 
tain circumstances. One useful approach to cope with 
non-Markovian situations is called the non-Markovian 
stochastic Schrodinger equation or the quantum-state- 
diffusion (QSD) equation [20, 21]; it describes an open 
system by a stochastic pure state rather than a density 
matrix. 

One of the successful cases that was dealt with using 
QSD is a system comprising two interacting qubits cou¬ 
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pled to a bosonic environment [22], which was studied by 
Zhao et al. They derived the exact time-local QSD equa¬ 
tion for such a system and discovered the entanglement 
generation caused by the environmental memory. How¬ 
ever, describing an open system with a reduced density 
matrix instead of a stochastic pure state would some¬ 
time be better since a stochastic state may carry an un¬ 
certainty represented by a complex stochastic Gaussian 
process. In contrast, the uncertainty in a reduced density 
matrix is ruled out by an ensemble average. In addition, 
in quantum control problems such as those discussed in 
[23], knowledge of the analytic form of the master equa¬ 
tion is quite important. Furthermore, since control prob¬ 
lems usually require a large amount of calculation, mak¬ 
ing a good approximation to reduce the computational 
complexity is a necessity. 

Based on the formal form of the master equation for 
multi-qubit dissipative systems [24], we have derived 
the analytic form of the master equation for this two- 
qubit system, where the noise is modeled by the Orn- 
stein-Uhlenbeck process. In particular, we write down 
an exact master equation composed of several functions 
whose time derivatives are analytically given. Moreover, 
we validate our analytic master equation by examining 
the entanglement generation and the purity of state, com¬ 
paring them to the results in [22]. In addition, we discuss 
the effects of first-order noise and show that from an 
asymptotic perspective, the exact master equation can 
be reduced to an approximate form, wherein the terms 
related to the first-order noise are eliminated. 

This study is organized as follows. In Sec. H, we intro¬ 
duce the model and the stochastic Schrodinger equation 
for the system comprising two interacting atoms coupled 
to a bosonic environment and derive the analytic form of 
the non-Markovian master equation. In Sec. HI, we use 
the non-Markovian master equation to examine the en¬ 
tanglement generation and state purity phenomena due 
to environmental memory and then validate the master 
equation. Next, we compare the solutions of the exact 
and approximate master equations in the last part of 
this section. Finally, Sec. IV presents a conclusion of 
this study. 
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II. MATHEMATICAL DESCRIPTION 
A. Model and stochastic Schrodinger equation 


For a system comprising two interacting two-level 
atoms coupled to a bosonic bath [22], the total Hamilto¬ 
nian can be divided into three parts (set fi = 1): 


/••• / ••• / (•), represents the en¬ 

semble average operation. 

In general, the functional derivative in the integral of 
Eq. (2) can be replaced with an operator 0{t, s, z*), such 
that 


SMz*) 

6zt 


0{t,s,z*)ii)t{z*). 


( 3 ) 


Hsys = + WsCrf -f Jxy{(T+(y- + 

+ erf, 

Hbath =yf 
i 

Hint = '^{gtalL + g*aiL^), 

i 


( 1 ) 


where wa and wb are the transition frequencies of 
two interacting atoms, a± = {ax ± iay)j2 is the cre¬ 
ation/annihilation operator for an atom, ai(aj) is the 
annihilation (creation) operator of the f**' mode of the 
bosonic bath, gt is the coupling constant between the sys¬ 
tem and the mode of the bath, and L = KAa{^ + KBa^ 
is a system operator. Note that the interaction between 
two atoms is modeled on the Heisenberg XXZ model, 
wherein the coupling constants Jx = Jy = Jxy are not 
necessarily equal to Jz- 

The non-Markovian stochastic Schrodinger equation 
for two interacting atoms coupled to a common bath is 
given as follows: 


For this two-interacting-qubit model, the O operator can 
be written as follows: 

0{t,s,z*) = Oo{t,s) + [ dsiz*^Oi{t,s,si), (4) 
Jo 


where Oo(t,s) and Oi(t,s,si) are related to the zeroth- 
and first-order noise components, which may be formu¬ 
lated as follows: 


Oo(t,s) 


Oi (t,s,si) 


fi{t, s)crf -f f 2 {t, s)a^ + h{t, s)a^a^ 

+ f4{t,s)aj:af, 
if5{t,s,si){2a^a^), 

( 5 ) 


in which the time derivatives and the initial conditions 
of the fj terms are listed in Appendix A. 


B. Exact non-Markovian master equation 

In general, the formal non-Markovian master equation 
can be derived from Eq. (2), and it reads as follows: 


^ipt =-iHsys'ipt +Lz^tp-L'^ j dsa{t,s)^, ( 2 ) 

where ipt = iptiz*) is a stochastic wave function of 
the system, and Zt = gjZje~^‘^^* is a complex 

Gaussian process satisfying A4[ztj = 0, A4[ztZsj = 
0, and A4[ztz*j = a(t,s), which is the correla¬ 

tion function of the bath and determines the envi¬ 
ronmental memory time. Here, the symbol, AI[*] = 


^=-t[H,y,,p,] + [L,M[PtO\t,z*)]] 

+ [M[0{t,z*)Pt],L^, 

where 0{t, z*) = dsa{t, s)0{t, s, z*), Pt = 

\tpt{z*)) {'ipt{z)\, and pt = M[Pt], which is the reduced 
density matrix used to describe the dynamics of the sys¬ 
tem. For a system that has the first order of noise as the 
highest order, A4[PtOf can be written as follows[24]: 


M[PtO^ = ptOlit) + [ dsi [ ds2a{si,S2)Oo{t,S2)ptO\{t,si) 
Jo Jo 


+ dsi dS 2 / dS 3 / dS4a{si,S2)a{s3,S4)*Oi{t,S2,S3)ptOl{t,S4)Ol{t,Si), 


( 7 ) 


where Oo{t) = J* dsa{t, s)Oo{t, s) and Oi(<, si) = 
Jp dsa(t, s)Oi (t, s, si). In our two-qubit system, the 
term involving the quadruple integral vanishes since 
OjOi = 0 due to cr+cr+ = 0. 

To make good use of Eq. (6), we have to convert Eq. (7) 
to a more explicit form, from which Ai[PtO^ can be 


easily computed. For this, we first define 


fjit) = [ ds a{t, s)fj{t, s), j = 1 - 4, (8) 

^0 
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and 

Fjit)= [ dsi / ds2 a{si,S2)fjit,S2)f5it,si), 

Jo Jo 

J = 1 - 4. 

When Oo and Oi given in Eq. (5) are substituted into 
Eq. (7), the equation gets converted to an explicit form 
as shown below: 

M[PtO^] = ) 

- 2i[Ficr:^ + F2a’^ + F^a'^a^ + F4a^af)pta^af. 

( 10 ) 

Here, we consider the Ornstein-Uhlenbeck process with 
the correlation function a{t,s) = which cor¬ 

responds to a Lorentzian power spectrum, to model the 
noise. In this process, with the help of the three other 
auxiliary functions 

Foit) = [ dsi [ ds2 a{si,S2)f5{t,S2)f^it,si), ( 11 ) 

^0 Jo 

hit, Si) = / ds a(t,s)/ 5 (t,s,si), (12) 

^0 

and 

hit)= f ds a{t,s)hit,s), 

Jo 


In fact, when the system is close to the Markovian 
regime, the operator Oi(t, s, Si) associated with the first- 
order noise in Eq. (4) becomes less important. To find the 
noise effect from Oi, we may approximate 0 (t, s, z*) with 
Oo(t, s) or equivalently let fsit, s, si) = 0 and then com¬ 
pare the results from the exact and approximate master 
equations. In the approximate master equation, M[PtO^ 
is reduced to 

M[PtO^]=Ptihai + f*af + f;a^af + f:ais^), (15) 

where the time derivatives of the fj terms form a set of 
coupled differential equations obtained by eliminating /s 
and /5 from Eq. (A5). 


III. NUMERICAL SIMULATIONS AND 
DISCUSSIONS 

In this section, we validate the newly derived master 
equation and compare the time evolutions of the density 
matrix pt governed by the exact and approximate O op¬ 
erators. Here, the states and operators are represented in 
the basis {| 11 ), | 10 ), | 01 ), | 00 )}, where CTz | 1 ) = | 1 ) and 
a. |0) = - |0). 


(13) A. Entanglement generation and state purity 


the explicit form of A4[PtO^ in Eq. (10) be¬ 
comes very useful, since and 

fo™ 1^0 sets of coupled dif¬ 
ferential equations (see Appendix A), implying that the 
analytic J^[PtO^ can be easily calculated by numerically 
solving these differential equations. With the derivation 
of Al[PtOt], we can compute pt from Eq. ( 6 ), and col¬ 
lectively, this completes the construction of the analytic 
form of the exact non-Markovian master equation for this 
two-interacting-atom system. 


C. Markovian regime and noise effects 

If the correlation function a(t, s) in Eq. (2) is given by 
7 i 5 (t — s), meaning that the quantum state at time t is 
independent of the state at time s (s < t), then there is no 
memory effect in the surrounding, and Eq. (2) becomes 

f) 

—■i/'t = -iHsysi’t + Lz^ipt - ^L^Lipt, (14) 

which is the stochastic Schrodinger equation in the 
Markovian limit. Note that in the Ornstein-Uhlenbeck 
process, when 7 —>■ 00 , the correlation function will be¬ 
have like a Dirac delta function, and hence, the environ¬ 
ment becomes memoryless. 


To validate the exact master equation, we examine 
whether it can reproduce the same entanglement gen¬ 
eration and time evolution of state purity as given by the 
QSD equation in [22]. Here purity is defined as tr[p^(t)], 
which conveys whether a system is in a pure or a mixed 
state. 

From [22], we know that if the initial state is i/'o = 110) 
(or po = | 10 ) ( 10 |), the two-qubit system evolves into a 
mixed and entangled steady state; while, if the initial 
state is V'o = :^(|11) + |00)), : 3 (| 10 ) + jOl)), or ]11), the 
system will first evolve to a mixed state because of the 
coupling interaction and finally decay to the pure state 
100 ) resulting from the dissipative environment. 

In Fig. 1, we plot the time evolution of purity for those 
four initial states generated by the exact master equation. 
We can see that unlike the other three initial states whose 
purities first become smaller than unity and then return 
to unity, the system remains in a mixed state with a 
purity value close to 0.5 only when ipo = |10). 

Moreover, from the time evolution of the density ma¬ 
trix with Po = 110) (10|, plotted in Fig. 2, we can clearly 
observe the system ending up in a state that is not only 
mixed but also entangled. Therefore, we can infer that 
the results from the exact master equation match well 
with the ones from the QSD equation described above, 
and this consistency conhrms the validity of our master 
equation. 




4 



FIG. 1. (Color online) Time evolution of purity for four differ¬ 
ent initial states. Here 7 = 1 , uja = ojb = 0.5, ka = bb = 1, 
Jxy = 0.7, and = 0.3. 



t 


t 


FIG. 2. (Color online) Dynamics of the density matrix with 
the initial state ipo = |10). Parameters are chosen as 7 = 1, 
OJA = ojb = 0.5, ka = Kb = 1, Jxy = 0.7, and Jz = 0.3. The 
values of other elements are always zero except pz 2 = p 23 - 


B. Exact master equation versus approximate 
master equation 


In the following, we assume that oja = ojb = <oJ and 
KA = Kb = K for simplicity. For such a condition, we can 
deduce from Eq. (A5) that ]\ = and /s = /4 given the 
symmetry in their time derivatives and the same initial 
conditions, share /i(0) = / 2 ( 0 ) = fsiO) = / 4 ( 0 ) = 0. 
Furthermore, if we substitute ijjt = Ci{t) |ll)-|-C2(f) |10)-l- 
czit) |01) -I- C 4 (t) |00) into Eq. (2), the dynamics of the Cj 


values are found to be 


—— — — i{ 2 ui Jz)ci — K( 2 fi + 2/3)01, 
df 

cic 2 ^ — — 

— = KZ*Ci - 2h:/5Ci -I- iJzC 2 - k(/i - / 3 )C 2 - iJxyCs 

- k(/i - /3)C3, 

cic3 — — 

— = KZ*Ci - 2K/5C1 -I- iJzCs - k(/i - / 3 )C 2 - iJxyC 2 

- k(/i - /3)C3, 


dC 4 

dt 


KZ*(c2 + C3) — i{Jz — 2 oj)c 4 , 


(16) 

where z*) = i Jq dsi/ 5 (t, si)z*^. Again, we can infer 
the relationship C 2 {t) = 03 ( 1 ) when C 2 ( 0 ) = 03 ( 0 ) because 
of the symmetry between C 2 and C 3 . For convenience, we 
shall assume C 2 ( 0 ) = 03 ( 0 ) without any loss of generality 
in the comparisons. 

Under these assumptions, pj 2 = M[cjC 2 ] = J^lcjC^] = 
Pj 3 (/ = 1 ^ 4), and thereby, together with the property 
Pt = pj of a density matrix, there are only six indepen¬ 
dent elements in pp pn, pi 2 , pi 4 , P 22 , P 24 , and P 44 . Next, 
we compare the differences of these six elements consid¬ 
ering the situations with or without the first-order noise. 

First, we consider the time evolution of a state, /'t, 
initialized as ijjo = |(|11) + |10) -I- |01) -I- |00)) with pa¬ 
rameters 7 = 0.1, oj = 0.5, and k = 1. From Fig. 3, we 
can infer that the dynamics of ijjt governed by the exact 
and approximate master equations only show very small 
differences although the environment is way off from the 
Markovian regime because of the small value of 7. Next, 
if we increase the value of 7 to unity, it is as expected 
that the system reaches its final state faster due to lim¬ 
ited back reaction from the environment in the Marko¬ 
vian regime. In this case, the results produced by both 
master equations match each other perfectly because the 
first-order noise is less important in the Markov limit. 
From Fig. 3, we can further deduce that the approximate 
master equation behaves well when oj increases from 0.5 
to 2, while the environment is still in a non-Markovian 
regime (7 = 0.1). In fact, variation of Jxy and Jz also 
does not initiate any clear difference between the results 
from both these master equations. 

Next, in Fig. 4, we show the dynamics from the same 
initial state governed by the exact and approximate mas¬ 
ter equations when 7 = 0.1, oj = 0.5, and k = 2. Com¬ 
pared with the red solid line and orange dashed line in 
Fig. 3 with K = 1, the oscillations are more intense with 
a shorter response time for each element of the density 
matrix, which can be explained by the analytical solution 
to Eq. (16). For example, from the solution for Ci{t) 


ci(t) = (17) 


we can infer that k has a direct influence on the oscilla¬ 
tions and the response time. Incidentally, the same effect 
can also be observed in the solutions to C2, C3, and C4. 
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FIG. 3. (Color online) Dynamics of the density matrix with different combinations of 7 and uj. Here, the initial state is 
V'o = |(|11) + |10) + |01) + |00))- Legend ( 7 , uj, e/a) indicates values of 7 and to for a curve, and e/a is used to indicate whether 
a curve is generated using the exact or approximate master equation. Other parameters are chosen as k = 1, J^y = 0.7, and 
Jz = 0.3. Note that the absolute value of each complex element is shown here. 



FIG. 4. (Color online) Dynamics of the density matrix with 
7 = 0.1, oj = 0.5, K = 2, Jxy = 0.7, and Jz = 0.3. The initial 
state is ■i/o = |(|11) + jlO) + |01) + |00)). 


Moreover, it is also observed that the results from the 
exact and approximate master equations do not match 
each other in the transient time. As for the cause, we first 
notice in Eq. (5) that the operator associated with the 
first-order noise is and then, apparently, the differ¬ 

ence comes from the coefficient ci{t) of state |11), which 
is operated on by a^a^. To verify this argument, we 
consider another initial state ipo = + |01) + 1^0)), 

where |11) is dropped out. With the same parameters 
that are used in the simulation shown in Fig. 4, we plot 
the dynamics of the new state without |11) in Fig. 5. It is 
evident from this figure that the results from both master 
equations as expected are well matched. 


Moreover, in Fig. 4, although the approximate dynam¬ 
ics do not mimic the exact ones in the transient time, they 
do match well in the steady state regime. Furthermore, 
even if k is increased to a higher value, the dynamics 
given by both master equations still match in the final 
state. 


From the above discussions, we can summarize the 
findings as follows: the master equation without the 
terms related to the first-order noise can be quite a good 
approximation, especially when state |11) does not form 
part of the initial state, or at least, the exact master 
equation can be asymptotically approximated by neglect¬ 
ing the effect from the first-order noise. Hence, we can 
also anticipate that residual entanglement and state im¬ 
purity from the initial state |10) would still be preserved 
without the first-order noise. 








































6 


0.4 r 


\P' 


'22 




■exact 

approx. 


0 10 20 30 40 50 60 70 


t 



1 


0.5 


\PiA 


0 -^^^^^-' 

0 10 20 30 40 50 60 70 

t 


IV. CONCLUSION 


In this study, we derived the analytic form of the non- 
Markovian master equation of a system comprising two 
interacting two-level qubits coupled to a bosonic envi¬ 
ronment. Following the derivation, the master equation 
was validated by examining entanglement generation and 
state purity. The exact and approximate master equa¬ 
tions were compared for various situations, and our re¬ 
sults reveal that the approximate master equation is well 
fitted to the exact master equation in steady state, and 
therefore, it serves as a good approximation in an asymp¬ 
totic sense. One of the potential applications of our re¬ 
sults is to implement two-qubit gates with this two-qubit 
system. Furthermore, the concept developed for this two- 
qubit system can be extended to multi-qubit systems. 
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Appendix A: Time derivatives of functions in the O operator and in the exact master equation 

The time derivatives of the coefficients fj (j = 1 ~ 5) for the O operator are as follows: 
dfi 

-^(t, s) = {2iuJA + KAfl + KBf3)fl + {-iJxy “ Ksfl + Kb/ 4)/3 + (2iJz + + Kb/ 3)/4 “ s), 


df2 

dt 

dh 

dt 

dfi 

dt 

dh 

dt 


{t,s) = {2iu}B + HaU + Kb/2)/2 + {-iJxy “ KAf2 + Ka/3)/4 + (2iJz + KaU + Kb/3)/3 “ il^Ah{t,s), 

{t,s) = {2iu}B + HaU + KBf2)f3 + {-iJxy “ Ka/2 + Ka/3)/i + (2iJz + HaU + «b/3)/2 “ iKAh{t,s), (Al) 
(t, s) = {2iu)A + KaIi + KBf3)f4 + {-iJxy “ Ks/l + KBfhh + {‘2.iJz + ^aJa + KbIs)!! “ il^Bhit, s), 

{t, S, Si) = [2i{uJA + Wb) + Ka(/i + fh + KB{f2 + f3)]h + («a/i + Ksh “ Kb/s “ KAfhfbit, Si), 


where fj{t) = /q ds a{t, s)fj{t, s) (j = 1 ~ 4) and / 5 (t,Si) = /g ds a{t, s)f 5 {t, s, Si). The initial conditions for these 
five functions are as follows: 

fi{t,t) = ka; hith) = i^b; f3it,t)=0; f4{t,t)=0; hit,'t,si) = 0; h{t,s,t) = -'i[KAf3{'t,s) + KBf4it,s)]. (A2) 
Next, we consider the time derivatives of fj{t) (j = 1 ^ 4), whose formal representations are expressed by 


dfj 

dt 


= a{t,t)fj{t,t) + J ds +y dsa(t,s) 


dfjit,s) 


dt 


= if 3 {t,t)-l dsa{t,s)fj{t,s)+ dsa{t,s) 


9hit,s) 

dt 


(A3) 


Note that the time derivatives_of / 5 (t,_si) and /^{t) have similar forms to Eq. (A3) with {fj{t, t),fj{t, s)} replaced by 
{hit,t,si),h{t,s,si)} and {hith), h{t, s)} respectively, where h{t,t) is given by 


h{th) = ds a{t,s)[-i{KAf 3 + Kb/ 4 )] = -i^KAh + Kb/ 4 )- 
do 

According to Eq. (A3), the time derivatives of fj (j = 1 ~ 5) and h are 

+ (-7 + 2za;^ + ka/i)/i + {-iJxy + kb/4)/3 + + kaU + kb/s)^ - iK.Bh{t, s), 

+ (-7 + 2iWB + KBf2)f2 + {-iJxy + Ka/s)^ + (2iJz + KaJa + KBf3)f3 “ iHA/bit, s), 


(A4) 


dt 

dt 

df5 

dt 

d /5 


(t) = (-7 + 2iuJB + KaJa + KBf2)f3 + {-iJxy “ Ka/2 + Ka/s)/! + (2zJz + KaJa + KBf3)f2 “ iKAhi^, «)) 

(t) = (-7 + 2iuJA + KaJi + KBf3)f4 + {-iJxy “ «:b/i + KBf4)f2 + (2iJz + KaJa + Kb/s)/! “ iKB/s)^, s), 

(t. Si) = [-7 + 2i{u3A + UJb) + 2KAfl + 2KBf2]h, 

.7, 


(A5) 


(i) = -i-^iKAh + kb/a) + [-27 + 2z(w^ + ujb) + 2KAfi + 2kb/ 2]/5- 
As for the time derivatives of Fj(t) (j = 1 ~ 4), we may also write down their formal representations as follows: 


dt 


= ds2a(t,S2)/i(t,S2)/5(t,t) + / dsia(si,t)/j(t,t)/5(t,si) 

^0 Jo 

+ dsi^ ds2a(si,S2)^(t,S2)/5(t,Si)dsi J ds2a(si, S2)/,(t, S2)^(t, Si), 


(A6) 


where the first two terms equate to fj{t)f^{t,t) and fj{t,t)fi{t) respectively. For dF^/dt, we just replace each fj 
term in Eq. (A6) with h- 



Below, we list the time derivatives of Fj{t) {j = 1 ~ 5): 


dFi 

dt 

dfz 

dt 

dfs 

dt 


dF4 

dt 

dfs 

dt 


iflinAfs + Kb ft) + KAft + (ka/i + KbH - 7 - 2za;B + 2KAft + 2kb/2)^’i + {-iJxy - kb/i + Kb/4)^3 

+ (2t+ KAf4 + KBh)Fi - IKbF^, 

ihiKAft + Kb ft) + Kb ft + (ka/4 + ^ 3/2 - 7 - 2tWA + 2KAft + ‘^KBft)F2 + {-iJxy - KAf2 + KAh)F4 

+ (2t+ KaJa + KBh)F3 - ikaFz, 

ih{KAft + Kb ft) + (kaIa + KBf2 - 7 - 2tWA + 2KAft + 2 ^ 5/2 )^3 + “ KAf2 + Ka/s)^! 

+ (2t+ KaJa + KBh)F2 - iKAFz, 

ihiKAft + Kb ft) + («a/i + KbH - 7 - 2tWB + 2.KAft + ‘2'KBft)FA + {-iJxy - KbJi + KBfA)F 2 
+ {2iJz + ka/4 + KBh)Fi - ikbF^, 

2Re[ih{KAft + «b/:)] + [-27 + 4KARe(/i) + 4KBRe(/2)]F5. 

(A7) 



